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Abstract. In 1963 Fisher and Stephenson [FS] conjectured that the monomer-monomer 
correlation on the square lattice is rotationally invariant. In this paper we prove a closely re- 
lated statement on the hexagonal lattice. Namely, we consider correlations of two quadromers 
(four-vertex subgraphs consisting of a monomer and its three neighbors) and show that they 
are rotationally invariant. 



1. Introduction 

The lines of the square grid divide the plane into unit squares called monomers. A 
dimer is the union of two monomers that share an edge, and a collection of disjoint dimers 
is said to be a dimer tiling of the planar region obtained as their union. We denote by 
M(i?) the number of dimer tilings of the region R. 

Let S n be the square [0, 2n] x [0, 2n]. Denote by S n (r, s;p, q) the region obtained from 
S n by removing the two monomers whose lower left corners have coordinates (r, s) and 
(r + p, s + q), < r, s, r + p, s + q < 2n — 1. 

The boundary-influenced correlation of the two removed monomers is defined in [FS] as 

M(S n (r,s;p,q)) 

u b {r,s;p,q) := hm — — . (1.1) 

If the monomers of S n are colored in a chessboard fashion, every dimer consists of 
one monomer of each color, and one sees that a>&(r, s;p,q) = unless the two removed 
monomers have opposite colors. Therefore in the following discussion we assume that the 
latter condition holds (this amounts to p and q having opposite parities). 

The correlation at the center is then defined as 

u(p,q):= lim u b (r,s;p,q). (1.2) 



r,s — >oo 
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The results in [FS] provide explicit exact values for u(p,q) for small, concrete values 
of p and q. In the cases q — (or q = 1) and p = q + 1 — which correspond to a lattice 
and a lattice diagonal direction, respectively — , Fisher and Stephenson [FS] computed the 
values of these correlations for the first thirteen (respectively, first 22) values of p. These 
tables provide strong evidence that the correlations u>(p,q) decay to zero at exactly the 
same rate (namely, as the inverse square root of the distance between the monomers, with 
exactly the same constant of proportionality) along this two inequivalent directions 1 . After 
making this remark, Fisher and Stephenson continue in [FS] by stating: "This equality 
suggests that the monomer correlations decay isotropically (with full angular symmetry)." 
This conjecture was the starting point of the present paper. 

Clearly, these considerations can also be made on the triangular lattice. A monomer 
is then a unit equilateral triangle, and the dimers are the unit rhombi consisting of the 
union of two monomers with a common edge. 

Each dimer covers one up-pointing and one down-pointing monomer, so the two orienta- 
tions of monomers play in this case the role of the two colors in the chessboard coloring of 
the square lattice. Therefore, the natural analog of the conjecture of Fisher and Stephen- 
son would be the statement that the correlation of two monomers of opposite orientation 
is rotationally invariant. 

A quadromer is by definition the union of any monomer with the three monomers 
adjacent to it. Clearly, this is just a lattice triangle of side two. 

It is correlations of two quadromers that we prove are rotationally invariant in this 
paper. Our method of proof is based on exact counting of dimer coverings of certain 
regions on the triangular lattice, and the reduction of the problem to these exact counts 
works only when we remove triangles of side two — quadromers — , and not when monomers 
are removed. In essence, this is due to the fact that two is an even number, and this allows 
to express the relevant quantities as convenient determinants, evaluate them and obtain 
the asymptotics of the resulting expressions. 

We note that removing a quadromer from a region is equivalent to removing any two 
of its outer monomers — such a pair could be called a birrier. Indeed, any dimer covering 
of the region resulting by removal of a bimer must contain the unit rhombus fitting in its 
notch, thus extending it to a quadromer. 

In order to define quadromer correlations, we need a family of regions on the triangular 
lattice to play the role the squares S n played on the square lattice. The choice of these 
regions is perhaps the single most crucial part in our proof. 

To this end, for any integer n > 1 we define the region P n as follows. Fix a vertex O 
of the triangular lattice and consider the zig-zag lattice path L u of length 2n extending 
upward from O and taking alternate steps northeast and northwest (see Figure 1.1). Let O' 
be the southwestern neighbor of O, and consider a second zig-zag lattice path of length 
2n, extending downward from O' and taking alternate steps southeast and southwest. The 
union of these two lattice paths and the segment 00' forms the eastern boundary of P n . 
The rest of its boundary is defined by moving from the lowest vertex of L d successively 



1 Hartwig [H] proved that the correlations decay as predicted by Fisher and Stephenson along the 
lattice diagonal direction. 
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Figure 1.1. P 4 . 



Figure 1.2. P 4 (3,0;2,l). 



n units west, 2n units northwest, 2n + 1 units northeast and n units east, along lattice 
lines throughout, to arrive at the topmost vertex of L u . Our region P n is the region traced 
out this way, with the additional requirement that the n dimer positions closest to L u 
have weight 1/2 (these positions are indicated by shaded ellipses in Figure 1.1): what this 
means is that a dimer covering of P n using precisely k of these dimers gets weight l/2 fe , 
and M(P n ) is the weighted count of the dimer coverings of P n , i.e., the sum of all their 
weights. 

Let i be the vertical line through O. Denote by D(R, v) the down-pointing quadromer 
whose base is centered R units to the left of £ and lies on a horizontal lattice line crossing 
Ld2v + 1 units below O' . Let U(R, v) be the up-pointing quadromer whose base is centered 
R units to the left of I and lies on a horizontal lattice line crossing L u 2v units above O 
(Figure 1.2 shows quadromers D(3,0) and £7(2, l)) 2 . 

Let P n (Ri, v±; R2, v 2 ) be the region obtained from P n by removing the quadromers 
D(R 1 , vi) and U(R 2 , 1*2), where R2 > 1 and v±, v 2 > (see Figure 1.2 for an example). 

Paralleling (1.1), we define the boundary-influenced correlation of two removed quadro- 
mers as 



In analogy to (1.2), we define the correlation of the quadromers at the center, with 



2 Considering the relative position to L u and L^, there are two distinct types of both up-pointing and 
down-pointing quadromers. The reason we made the indicated choice out of the total of four possibilities 
is that it provides even vertical separations between the bases of the quadromers. The other three choices 
would entail only minor alterations to the considerations in this paper. 



u)b{Ri, vi, R 2 , v 2 ) := lim 



M(P n (R u v 1 ;R 2 ,v 2 )) 
M(P n ) 



(1.3) 
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horizontal separation r and vertical separation \/3u > y/3, by 

uj(r,u):= lim Ub(R + r,u- 1;R, 0). (1.4) 

When making the separation (r, it) of the two quadromers grow to infinity it is natural 
to do it so that u = qr + c, where q and c are fixed rational numbers. 
The main theorem of this paper is the following. 

Theorem 1.1. As r and u approach infinity so that u = qr + c, with q > and c fixed 
rational numbers (c > 1 if q = 0), 

^)=4^VT3^ + ° (r_2) - (L5) 



Since the parenthesis in the denominator in (1.5) is just the square of the distance be- 
tween the centers of the bases of the removed quadromers, this theorem shows in particular 
that the quadromer correlation u(r,u) is rotationally invariant. 

Remark 1.2. An alternative (and perhaps more natural) way to define the correlation 
of the two removed quadromers would be to consider say regular hexagons H n of side n 
instead of our regions P n . One could even remove the two quadromers from near the center 
of H n: and then the analog of (1.3) would define directly the correlation at the center. 

However, when carrying out the line of approach of our proof in this set-up, instead of 
the well-behaved regions P n [ki, fo; h, fa] of Section 2, one is lead to regions obtained from 
H n by removing pairs of unit triangles from along two of its opposite sides — and these 
turn out not to possess a simple product formula for the number of their tilings. This 
precludes having an analog of Lemma 2.1, which is the first key step in our proof. 

2. A quadruple sum for quadromer correlations 

The reason we chose the regions P n as above is that they have, as shown in [C], the 
remarkable property that quite general alterations of their eastern boundary create regions 
whose weighted dimer counts are given by simple product formulas. 

To state this precisely, view the lattice paths L u and La on the eastern boundary as 
consisting of n bumps each — successions of two lattice steps forming an angle opening to 
the left (one sample bump on each of Ld and L u is indicated in Figure 1.1). Label the 
bumps on each path successively by 0, 1, . . . , n — 1, starting from the bumps closest to 
O 3 (see Figure 2.1). Let B be a bump on L d: and consider the unique down-pointing 
quadromer D that contains B. By removing bump B from P n we mean removing the 
three monomers of D contained in P n . We define removal of a bump of L u analogously, 



3 Note that this differs from the labeling used in [C], which is obtained by increasing the present labels 
by 1. 
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Figure 2.1. P 4 [0,2;l,2]. 



via the unique ^pointing quadromer containing it. Figure 2.1 illustrates the effect of 
removing bumps and 2 of and bumps 1 and 2 of L u from the region P4. 

Let P n [ki, k2',h, h] be the region obtained from P n by removing bumps k± and k 2 from 
L>d and bumps l\ and l 2 from L u , where < k± < k 2 < n — 1 and 0<Zi<Z 2 <n — 1 
(Figure 2.1 shows PJO, 2; 1, 2]). In [C, (2.2), (1.1)— (1.6)] explicit simple product formulas 
(i.e., with factors of size at most linear in the parameters) are given for the weighted count 
of dimer coverings of a family of regions R\ t q(x) that includes the P n [ki, k 2 ; h, foY s (1 and 
q are lists of strictly increasing positive integers, and x is a nonnegative integer; in the 
notation of [C], P[i,...,fc 1 ,fc 1 +2,...,fc 2 ,fc 2 +2,...,n],[i,...,/ 1 ,/ 1 +2,...,/ 2 ,/ 2 +2,...,n]H is the region that we 
denote here P n [ki, k 2 ; h, l 2 ]). Thus we obtain simple product formulas for the numbers 
M(P n [ki, k 2l /1, l 2 }). Using them, one obtains after straightforward if somewhat lengthy 
manipulations that 



lim 

n^oo 



M(P TC [fci, fc 2 ; IjM) ^ 2- 4 (2h + 1)! (2k 2 + 1)! (2h + 1)! (2/ 2 + 1)! 
M(P n ) 2 2 ^+ 2k ^+ 2 ^+ 2 ^k 1 l (h + 1)1 k 2 \ (k 2 + 1)1 hi 2 1 2 \ 2 

x (fc 2 -fci)(/ 2 -/i) 

X (h + h+ 2)(k t + l 2 + 2)(k 2 + h + 2)(k 2 + l 2 + 2)' ^ ^ 
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Lemma 2.1. The boundary-influenced correlation Ub(Ri, i>i; R 2 , v 2 ) is given by 



u h {R u Vl - R 2 , v 2 ) =2- 4 R t R 2 (R 2 - 1/2)(R 2 + 1/2) 

.1 ST, ua+b+c+d (fli + a-l)!(fli + &-!)! 





(R 2 + c-iy.(R 2 + d-iy. 



(2c+l)\(R 2 -c)\(2d+l)\(R 2 -d)\ 

(2ui + 2a+ 1)! (2ui + 26+l)! 



2 2(2« 1 +a+6)( t , l + a )| ( Vl + a + 1)1 ( Vl + ft)! + 6 + 1)! 



(2i> 2 + 2c + l)!(2i> 2 + 2d+l)! 

22(2«2+c+d)( V2 + c )!2( V2 + rf )|2 



(6-a) 2 (rf-c) 2 



(2.2) 



(u + a + c)(u + a + d)(u + b + c)(u + b + d) 



where u = v\ + v 2 + 2. 

To prove this Lemma we will need the following special case of the Lindstrom-Gessel- 
Viennot theorem on non-intersecting lattice paths. 

Our lattice paths will be paths on the directed grid graph Z 2 , with edges oriented so 
that they point in the positive direction. We allow the edges of Z 2 to be weighted, and 
define the weight of a lattice path to be the product of the weights on its steps. The weight 
of an iV-tuple of lattice paths is the product of the individual weights of its members. The 
weighted count of a set of iV-tuples of lattice paths is the sum of the weights of its elements. 

Let u = (u±, . . . , un) and v = (v±, . . . ,vn) be two fixed sets of starting and ending 
points on Z 2 , and let A/"(u, v) be the set of non-intersecting lattice paths with these 
starting and ending points. For P e A/"(u, v), let op be the permutation induced by P on 
the set consisting of the iV indices of its starting and ending points. 

Theorem 2.2 (Lindstrom-Gessel-Viennot [GV]). 



where aij is the weighted count of the lattice paths from Ui to Vj. 

What makes possible the use of this result in our setting is a well-known procedure 
of encoding dimer coverings by families of non-intersecting "paths of dimers:" given a 
dimer covering T of a region R on the triangular lattice and a lattice line direction d, the 
dimers of T parallel to d (i.e., having two sides parallel to d) can naturally be grouped into 
non-intersecting paths joining the lattice segments on the boundary of R that are parallel 
to d, and conversely this family of paths determines the dimer covering (see Figure 2.2 for 
an illustration of this and e.g. [C] for a more detailed account). 

We will find it convenient to view the paths of dimers directly as lattice paths on Z 2 , 
thus bypassing the "extra steps" of bijecting them with lattice paths on a lattice of rhombi 



(-l)' Jp wt(P)=det((a u ) 



l<i,j<n) -i 



PeAf(u,v) 
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Figure 2.2. Figure 2.3. 

A dimer covering of -P4(3, 0; 2, 1) The effect of Laplace expansion over 

encoded by paths of rhombi. the rows indexed by a and (5. 



with angles of 60 and 120 degrees, and then deforming this to the square lattice. In this 
context, the "points" of our lattice C are the edges of the triangular lattice parallel to a 
chosen lattice line direction d — we call them segments — , and the "lines" of C are sequences 
of adjacent dimers extending along the two lattice line directions different from d. The 
dimers that these lines consist of are the edges of C. 

Proof of Lemma 2.1. Choose the lattice line direction d in the above encoding procedure 
to be the southwest-northeast direction, and choose the positive directions in the lattice C 
so that they point east and southeast. Encode the dimer coverings of P n (Ri, vi; R2, v^) by 
(2n + 3)-tuples of non-intersecting paths consisting of dimers parallel to d (see Figure 2.2; 
there and in the following figures the dimer positions weighted by 1/2 are not distinguished, 
but are understood to carry that weight). 

Let P be such a (2n + 3)-tuple. Consider the permutation a-p induced by P on the 
set of the 2n + 3 indices of its starting and ending points. We claim that the sign of op 
is independent of P. Indeed, denote by U and D the removed up-pointing and down- 
pointing quadromers, respectively. While the way in which the starting and ending points 
of P — clearly independent of P — are matched up depends on P, it is always the case that 
the two paths ending on U start at consecutive starting points, and the two starting at D 
end at consecutive ending points. It is easy to see that this implies that all cr P 's have the 
same sign. 

Weight by 1/2 the edges of our "path-encoding" lattice C corresponding to dimer po- 
sitions weighted by 1/2 in P n (i?i, v\\ R2, V2). Weight all other edges of £ by 1. Then the 
weight of the dimer covering encoded by P is just wt(P), and we obtain by Theorem 2.2 
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and the constancy of the sign of op that 



M(P n (R u v 1 ;R 2 ,v 2 )) = | det A | 



(2.3) 



where A is the (2n + 3) x (2n + 3) matrix recording the weighted counts of the lattice paths 
with given starting and ending points (note that the right hand side of (2.3) is independent 
of the ordering of these starting and ending points). 

We deduce (2.2) by applying Laplace expansion to the determinant in (2.3). Recall 
that for any mx m matrix M and any s-subset S of [m] := {1, . . . , m}, Laplace expansion 
along the rows with indices in S states that 



det M = J2(-iy (K) det Mg det Mf™ 1 ^, 



(2.4) 



K 



where K ranges over all s-subsets of [m], e(K) := J2 keK (k — 1) and Mj is the submatrix 
of M with row-index set / and column-index set J. 

The rows and columns of the matrix A in (2.3) are indexed by the starting and ending 
points of the (2n + 3)-tuples of non-intersecting lattice paths encoding the dimer cover- 
ings of P n (Ri, v i; i?2, v 2 ). The starting points are the 2n + 1 unit segments along the 
northwestern boundary of P n (Ri, v±; R2, V2), together with the two segments a and (3 of 
D parallel to d (see Figure 2.2). The ending points are the 2n + 1 segments parallel to d 
on the eastern boundary of P n (R\, vi; R2, V2), together with two more such segments on 
U. 

Apply Laplace expansion to the matrix A of (2.3) along the two rows indexed by a and 
(3 (see Figure 2.2). The first determinant in the summand in (2.4) is then just a two by 
two determinant. Its entries are weighted counts of lattice paths on C that start at a or 
P and end at some segment on L d . There are only Ri + 1 segments on L d that can be 
reached this way. Label them consecutively from top to bottom by 0, 1, . . . , Ri. We can 
restrict summation in (2.4) to the two-element subsets K of this set of segments: all other 
terms have at least one zero column in the two by two determinant. Therefore we obtain 
from (2.3) that 



M(P n (R 1 ,v 1 ;R 2 ,V2)) 



£ (-l) a+& detA{:5 } detA|-^j\{:5 ; 



0<a,b<iti 



(2.5) 



Choosing the origin of C to be at a, one sees that f3 has coordinates (—1,1) and the 
segment labeled j on L d has coordinates (Ri — 1 — j, 2j), j = 0, . . . , R±. Since the lattice 
paths counted by the entries of ^j^,'^ have all steps weighted by 1, the determinant of 
this matrix is 



det A 



{a,b} 
{a,f3} 



r (Ri 



det 



(R!-l+a\ (R x -l+b\ 
\ 2a A 26 ) 

1 A 21 



2a- 



26-1 



2#! 

8 



(fr-a)QRi + a- 1)1 (Ri + b- 1)! 
(2a)!(i2i-a)!(26)!(i?i-6)! ' 



(2.6) 



On the other hand, the second determinant in the summand in (2.5) can be interpreted 
as being the weighted count of dimer coverings of the region P^ ,b \R 2 ^v 2 ) obtained from 
P n (Ri, vi; R 2 , 1^2) by placing back quadromer D and removing the two monomers (i.e., 
unit triangles) that contain the segments labeled a and b on L^, under the labeling of the 
preceding paragraph (see Figure 2.3 for an illustration). Indeed, the Lindstrom-Gessel- 
Viennot matrix of this region is precisely ^[2n+3]\{a 0}' anc ^ the argument that proved 

(2.3) we obtain that M(Pn ,b ^ (R 2 , v 2 )) is equal to det ^[2n+3]\{a py U P ^° a s *S n ^at is 
independent of a and b (indeed, the permutations op that occur when applying Theorem 
2.2 to the region pj\, ' ^ (-R2? ^2) f^re independent of a and 6). Therefore, using (2.6) we can 
rewrite (2.5) as 



M{P n {R 1 ,v 1 ;R 2 ,v 2 )) 



2R X 



V (-ir+ b {b ~ a){Rl + a ~ 1)! {Rl + b ~ 1)! M(P^(R 2 v 2 )) 
o<-tW (2a)l(R 1 -a)l(2b)l(R 1 -b)l ( - 



(2.7) 



In turn, M(Pt M (R2, v 2 )) can be expressed by a formula similar to the one above. To 
obtain this, encode the tilings of Pn' b \R 2 ,v 2 ) by lattice paths, choosing this time the 
lattice direction d to be the southeast-northwest direction, and the positive directions in 
the encoding lattice C to point east and northeast. As in the previous encoding, weight by 
1/2 those segments of C that correspond to dimer positions weighted 1/2 in p}f l ' b \R 2 , v 2 ), 
and weight all its remaining segments by 1. 

Each tiling of P^ 1,b \R 2l v 2 ) gets encoded this way by a (2n + 2)-tuple P of non- 
intersecting lattice paths on C, starting at the unit segments on its southwestern boundary 
or at the unit segments 7 and 6 of U that are parallel to d, and ending at the unit seg- 
ments parallel to d on its eastern boundary (see Figure 2.4). As in the argument that 
proved (2.3), the sign of the permutation op is independent of P. Therefore, we obtain 
by Theorem 2.2 that 

M(Pi a ^(R 2 ,v 2 )) = edetB, (2.8) 

where B is the (2n + 2) x (2n + 2) matrix recording the weighted counts of the lattice 
paths with specified starting and ending points, and the sign e in front of the determinant 
is the same for all choices of a and b. 

Apply Laplace expansion in det B along the two rows indexed by 7 and 5. The first 
determinant in the summand of (2.4) is again two by two, and records weighted counts of 
lattice paths starting at 7 or 5 and ending at some segment on L u (see Figure 2.4). There 
are R 2 + 1 segments on L u that can be reached this way; label them consecutively from 
bottom to top by 0, 1, . . . , R 2 . As with our previous Laplace expansion, we can restrict 
the summation range in (2.4) to obtain 



M{Pt b \R2,v 2 )) = e J2 (-l) c+rf - X det det ^g^^^;^- (2.9) 

0<c<d<R 1 



Figure 2.4. Figure 2.5. 

[1 21 

A dimer covering; of PI ' J (2, 1) The effect of Laplace expansion over 

encoded by paths of rhombi. the rows indexed by 7 and S. 

Centering C at 7, 5 has coordinates (—1, 1) and the segment labeled j on L u has 
coordinates (R 2 — 1 — j, 2j + 1), j = 0, 1, ... , R 2 . The weighted counts involved in the 
entries of B^ d ^ (which involve this time some steps of weight 1/2) are easily calculated 
and one obtains 



detS { { ;;f } =det 



l(R 2 -l+c\ , (R 2 -l+c\ l(R 2 -l+d\ , (R 2 -l+d\ 

2\ 2c ) ' \ 2c+l ) 2\ 2d ) ' \ 2d+l I 

l(R 2 -l+c\ , (R 2 -l+c\ l(R 2 -l+d\ , (R 2 -l+d\ 

2\ 2c-l J f l 2c / 2^ 2d-l / V 2d ) 

(d-c)(R 2 + c-l)\(R 2 + d-l)\ 



--2R 2 (R 2 -l/2)(R 2 + l/2) 



(2c + 1)! (R 2 - c)\ (2d + l)\(R 2 -d)\ (210) 



On the other hand, by applying Theorem 2.2 one more time one sees that 

det5[ 2 2 :^ (2.11) 

where pt h][cA is the region obtained from pjf'^lr^, v 2 ) by placing back quadromer U 
and removing the two monomers near L u containing segments c and d, and the sign e' is 
independent of c and d. Furthermore, p^' b ^ c,d ^ differs from the region P n [vi+a, vi+b; v 2 + 
c, v 2 + d] considered at the beginning of this section only in that the former contains four 
more dimers, which are weighted by 1 and forced to be part of all its dimer coverings (see 
Figure 2.6); so the two regions have equal weighted counts of dimer coverings. Therefore, 
by (2.7), (2.9), (2.10) and (2.11) we obtain that 
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Figure 2.6. 



The regions p[ h2][0 > 1] and P 4 [l, 2; 1, 2] 



differ only by four dimers. 



M{P n {R 1 ,v 1 ;R 2 ,v 2 )) = AR X R 2 {R 2 - 1/2) (R 2 + 1/2) 





(d-c)(fl2 + c-l)!(fl 2 + d-l)! 



(2c+l)!( J R 2 -c)!(2d+l)!( J R 2 -rf)! 
x M(P n [i;i + a,ui + b; v 2 + c, v 2 + d})\. 



(2.12) 



Dividing (2.12) by M(P n ), letting n — > oo and using (2.1), one obtains an expression for 
Ub(Ri, vi;R 2 ,v 2 ) as a quadruple sum in which the summation indices need to satisfy a < b 
and c < d. 

The fortunate situation is that, on the one hand, when a = b or c = d the summand 
becomes zero, and on the other, (2.12) and (2.1) combine to produce a summand which is 
invariant under independently transposing a with b and c with d (because the differences 
b — a and d — c end up appearing at the second power). Therefore the summation range 
may be extended to the one shown in (2.2), at the expense of a multiplicative factor of 
1/4. This leads precisely to the quadruple sum given in the statement of the Lemma. □ 



A simple partial fraction decomposition of part of the summand in (2.2) affords a great 
deal of simplification in the expression (2.2) of the boundary-influenced correlation. 



3. Four double sums, with integral representations 
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Indeed, one readily checks that 



(b — a)(d — c) 



(it + a + c){u + a + d)(u + b + c){u + b + d) 

1 



(u + a + c)(u + b + d) (u + a + d)(u + b + c) ' 

Using this, the portion of the summand contained in the last line of (2.2) becomes 

(b-a)(d-c) (b-a)(d-c) 
{u + a + c){u + b + d) (u + a + d)(u + b + c) 

ac + bd — ad — be ac + bd — ad — be 



(u + a + c)(u + b + d) (u + a + d)(u + b + c) 



(3.1) 



Using this, the fourfold sum of (2.2) becomes a difference of two fourfold sums. The 
advantage of this expression is that each of the latter two fourfold sums can be written 
as the product of two double sums. Indeed, in all the factors of the summand in (2.2) 
except the last line of (2.2), the summation variables can be separated. Furthermore, the 
variables {a, c} can be separated from {6, d} in the first term on the right hand side of 
(3.1), while {a, d} can be separated from {6, c} in the second term of the right hand side 
of (3.1). The double sums arising this way are all of the form 



Ri R 2 

a=0 c J~^ a C (2a)! (#i - a)! (2c + 1)! (R 2 - c)\ 

2t>i + 2a + l)! 2i! 2 + 2c + l)! 



M„(R U R 2 ).-^JJ-1) {2a) u Ri _ a) u 2c+1)l{R2 _ c)l 



x 



where v has one of the values 1, a, c and ac. For notational convenience we will often 
write simply M v instead of M v (Ri, R 2 ). 

More precisely, consider the term ac of the numerator of the first fraction on the right 
hand side of (3.1). When summing over a, 6, c and d as required by (2.2), this term gives 
rise to M ac M\. Similarly, the remaining terms of that numerator, bd, —ad and —6c, give 
rise to M\M ac , —M a M c and — M a M c , respectively. In the same fashion, the second term 
on the right hand side of (3.1) generates the products —M a M c , —M a M c , M ac M\ and 
MiM ac , respectively. Therefore, we obtain by (2.2) that 

u b (R u v i; R 2 , v 2 ) = 2~ 4 R 1 R 2 (R 2 - 1/2){R 2 + l/2)|4MiM ac - AM a M c \ 

= 2~ 2 R 1 R 2 {R 2 - 1/2)(R 2 + l/2)|MiM ac - M a M c \. (3.3) 

Thus, the asymptotic study of the correlation reduces to studying the asymptotics of these 
four double sums. By their definition (3.2), had it not been for the factor l/(u + a + c), 
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these double sums would further be separable as products of simple sums. We can get 
around the obstacle posed by this factor by expressing it as an integral 4 : 



1 



= / x 



u+a+c— 1 



dx. 



u + a + c 

Indeed, substituting this into (3.2), we obtain that for instance Mi is expressed as 

Ri 



(3.4) 



M MVV (fli + o-l)l (2^ 1 + 2a + l)! 



2 2 ^+ 2a (v 1 + a)\ (vt + a + iy. 
(R 2 + c-l)\ (2v 2 + 2c+l)l 



x" 



v c=0 



(2c + 1)! (R 2 - c)\ 2 2 ^+ 2c (v 2 + c)! 2 



x c x u 1 dx. 



(3.5) 



It is not difficult (indeed, with access to a computer algebra package like Maple, it is 
immediate) to see that the two sums in the above integral can be expressed in terms of 
hypergeometric functions 5 as 



y\ iy » (Ri + a-l) 



{2vx +2a + 1)! 



a=0 



(2a)! - a)! 2 2 ^+ 2a ( Vl + a)! ( Vl + a + 1)! 



J_ (2^ + 1)! 
i?i 2 2 ^i! (ui + 1)! 



5-^2 



i^i + 2 5 4 



(3.6) 



and 



R2 



E(-D' 



c=0 



QR 2 + c-l)! (2i> 2 + 2c+l)! 
(2c + l)!(i22-c)! 2 2 ^+2c( Vl +c )!2 

1 (2u 2 + l)! 



X ' 



R 2 2 2v *v 2 \ 2 



3 r 2 



—R 2 , R 2 , f 2 + o x 



,v 2 + l 



(3.7) 



Replacing these formulas in (3.5) we obtain the following result. 



4 This useful trick was pointed out to the author independently by Ira Gessel and Doron Zeilberger. 
Its use allows a shorter proof for the asymptotics of the four double sums than our original proof, which 
relied on the expansion 



J_ = I/i_i±£ + 

a + b u { u 



(a + c)(a + c-l) (a + c)(a + c- l)(a + c-2) 
u + a + b u I u u(u + 1) u(u + l)(u + 2) 

5 The hypergeometric function of parameters a±,.. . ,a p and b\, . . . , b q is defined by 

ai, ■ ■ ■ , a p 



+ ■ 



}• 



pF q 



bi, ... ,6c 



(ai)fc • ■ ■ japh z k 

^ k\(b!) k -.-(b q ) k ' 



where (a)rj := 1 and (a)*, := a(a + 1) • • • (a + k — 1) for fc > 1. 
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Proposition 3.1. The double sum Mi has the integral representation 
1 (2ui + 1)! (2u 2 + 1)! 



Mi 



R X R 2 2 2 ^+ 2 ^v 1 l(v 1 + l)lv 2 \ 2 



I :\T> 
'0 



-Rl, R\, V\ + s •<' 



3 



3-^2 



— i? 2 , -^2, ^2 
3 
2 



§,V2 + 1 



2 . 



X 



' 4 



(3.8) 



When applying the same reasoning to the remaining double sums M a , M c and M ac , 
two more sums besides (3.6) and (3.7) need to be expressed in terms of hypergeometric 
functions. These are 

"' in (fli+a-l)! 



a=0 



2 Vl + 2a + l)! 
(2a)! - a)! 2 2u i+ 2a (f 1 + a)! + a + 1)! 



x a 



(2t;i + 3)(2t; 1 + l)l 
1 2 2 ^+V!(^i + 2)! 3 2 



+ + l,Ui + 

§,vi + 3 



! 4 



(3.9) 



and 



i?2 



c=0 



CR2 + C-I)! 2i> 2 + 2c + l)! 
(2c + 1)! (R 2 - c)\ 2 2 ^+ 2c ( Vl + c)! 2 



x c c 



R 2 x{2v 2 + 3)(2v 2 + iy. 



3-^2 



-i? 2 + 1,R 2 + 1,V2 + 



,^2 + 2 



(3.10) 



3 2^+2^1 ( V2 + 1)1 

Using (3.2), (3.4), (3.6), (3.7), (3.9) and (3.10), we obtain the following result. 
Proposition 3.2. The double sums M a , M c and M ac have the integral representations 



M n 



R x (2i> 1 + 3)(2i> 1 + l)!(2i> 2 + l)! 
R 2 2 2 ^+ 2 ^+ 2 v 1 \(v 1 + 2)\v 2 \ 2 



x / .•',/• 2 
'0 



-R x + + l,vi + f _ x 



|,vi + 3 



3-^2 



-i? 2 , -R2, v 2 + § _ a; 
|,V2 + 1 ; 4 



(3.11) 



i? 2 (2i> 1 + l)!(2^ 2 + l)!(2i'2+3) 
" 3i2l 22^+2^+2^1 ( Vl + x)! v 2 ! (v 2 + 1)! 



x 



Jo 



3F2 



-Ri, Ri, v\ + 2 
i^i + 2 ; 4 



3-^2 



#2 + 1,^2 + 1,U2 + § . X 



l,v 2 + 2 



'4 



(3.12) 



R X R 2 (2 Vl + 3)(2m + 1)! (2i> 2 + 3)(2t> 2 + 1)! 
3 2 2 ^+ 2 ^+ 4 v 1 \(v 1 + 2)\v 2 l(v 2 + iy. 
.1 



/ sF 2 
'0 



-^1 + 1,^1 + 1,^1 + 1 
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|,vi + 3 



3-^2 



-^2 + 1,^2 + 1,^2 + 1 a; 



1,^2 + 2 



' 4 



(3.13) 



4. Jacobi polynomials 

The siVs of the preceding section can be expressed in terms of 2-Pi's by the formula 



3-^2 



a, 6, c 
a — n,d 



(i 



(b)k(c) k k j? 
(d)k 



b + k, c + k 
d + k 



^ (4-1) 

where n is a nonnegative integer. This follows for instance by comparing the coefficients 
of the powers of z and applying the Chu-Vandermonde summation (for the latter, see e.g. 
[GR, (1.2.9), p.2]). 

Applying this to the 3-F 2 's in (3.6) and (3.7), one obtains 



3-^2 



--Ri,-Ri,vi + 



x 



|,^i + 2 '4 



2~M + 1 



k=o 



k J ( "^"2 k+1 - fc (t, 1+ 2) fc 4^ 2Fl 



vi + 2 + fe '4 



(4.2) 



and 



3^2 



— -R27 -R27 ^ 2 + 9 ^ 



|,^2 + 1 '4 



(-V2 ~ \) 



V2 



X 



£(-u 



-#2 + l,R 2 + l.x 
V2 + I + I '4 



In turn, the resulting 2 Fi's can be expressed in terms of the Jacobi polynomials P, 
using the formula (found for instance in [Sz2, (4.21.2), p. 62]) 



2-^1 



n + a + /? + 1, —n _ 1 - x 
1 + a ' 2 



nir(l + a) (at0) 



r( 



(4.3) 

(a '«(x) 

(4.4) 



We obtain for the 2 Fi's of (4.2) and (4.3) the expressions 



-Ri + k, Ri + k _ x 
V! + 2 + k '4 



(gi - fc)! (m + fc + 1)! (vi+fc+ i, fc _ t , 1 _ 2 ) / x\ 

(i? 1 + ^ + l)! I 2y> (4 ' 5j 



and 



-R 2 + l,R 2 + l . £ 

t; 2 + 1 + / '4 



(i?2 -/)!(t>2 + /)! („ 2+ /,/-„ 2 -l) / _ X\ , . 

(^2 + ^2)! I 2)- {4 ' b) 



Substituting (4.2), (4.3), (4.5) and (4.6) in the integral representation (3.8) of M 1; we 
obtain that 



Mi 



Ul + l V 2 

EE 



CfcZ (-Ri)k(Ri)k 



{-R2)l{R2)l 



j^f Q J—^ R1R2 (Rl — k + l)„ 1+ fc+i (i?2 — I + l)v 2 +l 
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Iki{RiiR2)i 



(4.7) 



where Cki depends only on k, I, v\ and V2, and for < k < vi + 1 and < / < V2, 







p ( Vl +k+l,k- yi-2) 



(i-^tr-'H'-S 



dx. (4.8) 



By (1.4) and (3.3), we need the asymptotics of M\ for v\ = u — 1, = 0, R\ = R + r and 
i?2 = -R, where r > and u > 1 are fixed and — > oo. 

By (4.7), the study of this asymptotics of Mi reduces to the asymptotics of the /fe/'s. 
The asymptotics of the Jacobi polynomials pjf'^ for large n is given by the Darboux 
formula (see e.g. [Sz2, Theorem 8.21.8, p.196]) 



a + (3 + 1 



cos 



D^\cos9) = 



n + 



9- 



nil \ sin - 



a 




2 




e\ 




2) 





TV 



(4.9) 



E^\cos9) =0(n~ 3 / 2 ), 

where «,/3eR and < 9 < tt are fixed. 

We show in Lemma 4.2 below that replacing in I^i the Jacobi polynomials by their 
Darboux approximations leaves the asymptotics of Iki unchanged. This will follow from 
the following inequalities due to Szego [Szl] (see also [Sz2, p. 197]). 

Proposition 4.1 (Szego [Szl, (46'), (48'), p.77] [Sz2]). (a) Let a > -1/2, (3 6l and 
e > 0. Then there exists a constant A depending only on a, (3 and e so that 



P^\x) 



< 



A 



1 



\/n(\- x)f +3 



■1 + e < x < 1. 



(4.10) 



(b) Let a, (3 G R and e, c > 0. Then there exists a constant B depending only on a, (3, e 
and c so that 

B 1 



E^\cosd) 



< 



n 3/2 r +| ' 



cn 1 < 9 < 7r — e. 



One readily checks that 



cos _1 (x) > v^Vl -x, < x < 1. 



(4.11) 



(4.12) 



Indeed, if/(x) := cos" 1 ^) - v / 2v / T T7 ^, one has /'(x) = (l-x)- 1 / 2 (2" 1 / 2 -(l+x)- 1 / 2 ) < 
0, for x G [0, 1), and as jf(l) = 0, (4.12) follows. 

Using (4.12) and Proposition 4.1(b) one obtains that for any e > there exists constants 
B e and N e depending only on a, (3 and e so that 



E^\x) 



< 



B f 



1 



n J/ (1 - x) 2- + 4 2 
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(4.13) 



By Proposition 4.1(a), there exists a constant A depending only on a and (3 so that 



< 



A 



(i- x y 



< X < 1. 



(4.14) 



Using (4.9) and the fact that sin(cos _1 (x)/2) = a/(1 — x)/2 and cos(cos _1 (x)/2) = 
a/(1 + x)/2, one obtains that the Darboux approximants satisfy an inequality of the same 
form as (4.14): there exists a constant A' such that 



D^\x) 



< 



A' 



Vn (1 - £)f+i 



- < X < 1. 

2 - ~ 



(4.15) 



By (4.14) and (4.15), E^ ,P) {x) = Pt' P \x) - D { n a ^\x) also satisfies an inequality of the 
same type, so there exists a constant A" depending only on a and f3 so that 



X 



< 



A' 



(i-xy 



< X < 1. 



(4.16) 



We are now ready to prove the announced invariance of the asymptotics of Iki(Ri, R2) 
under replacement of the Jacobi polynomials by their Darboux approximants. Define 

JuiRuRi) := £ D %+W-<*-V (l _ |) x>(£M.'-<»-i) _ |) ^-H.-^. (4 . 17) 

Lemma 4.2. Let r, v\, V2, < k < v\ + 1 and < / < V2 be fixed, and assume that not 
all ofv\, V2, k and I are zero. Then we have 



lim R (I kl (R + r,R)- J kl {R + r, R)) = 0. 

R— >oo 



Proof. Denoting for simplicity f(x) := P^^ k k+1,k Vl (x), g(x) := Pr^-i' 1 V2 
F(x) := D { ^ +1 > k - Vl - 2 \x) and G(x) := D£+ M " t,2_1) (a;), we obtain 

\Iki{Ri, R2) — JkiyRii -R2) I 

= / {/(l - x/2)g(l - x/2) - F(l - x/2)G(l - x/2)} x^+^dx 

= I / {/(l - x/2)(/(l - x/2) - /(l - x/2)G{l - x/2) + f(l - x/2)G(l - x/2) 
Jo 

- F(l - x/2)G(l - x/2)}x k+l+u ~ l dx\ 

< [ |/(l-x/2)| \g(l- x/2) -Gil-x^lx^+^dx 
Jo 

+ [ \G(1 - x/2) I |/(1 - x/2) - F(l - x/2)G(l - x/2)\x k+l+u ~ 1 dx. 



Therefore, to prove the Lemma it suffices to show that for Ri = R + r and R 2 = R, 



lim R f pg+t +1 ' k -*-*\l-x/2) 
•R^ 00 Jo 



E 



(v 2 +l,l-v 2 -l) 

R 2 -l 



(l-x/2) x^+^dx = (4.18) 



and 



lim Rf E^ k k +1 ^- 2 \l-x/2) 
R^°° Jo 



< 2 i M """ 1) (i-^/2) 



x k+l +u - ldx = 

(4.19) 



By the Cauchy-Schwarz inequality, for h with integrable square we have 

h(x)dx^j < J h 2 (x)dx. 
Let < e < 1 be arbitrary. By (4.20) we have 



(4.20) 



(/V£-^ +1,fc """ 2)(1 - x/2) 



E 



< 



1-e 



(v 2 +l,l — v 2 — 1) 
R2-I 

2 



(l-x/2) 



X 



k+l+u-1 



dx 



p ^ 1 +fc +1)fc _ t)1 _2) (1 _ ^ ^+M-«2-l) (1 _ x/2) 



Jl-e 



R 2 - 



j^{v 2 +l,l-v 2 -l) 

R 2 —l 



X 



2(k+l+u-l) ^ x 



(l-x/2) 



x ^(k+l+u-l)^ x _ 



(4.21) 



Denote the two integrals on the right hand side of (4.21) by I\ and I2, respectively. By 
(4.13) and (4.14), we have 

-1-e 



A 2 B 2 



h < 



(R 1 -k)(R 2 -ir Jo 



{x/2)-^ +k+1+ ^{x/2)-^ +l+ ^x 2 ^ l+u -^dx 



< Me f 1 ^ T Vl+V 2 +k+l-li 

" RiRl Jo 



(4.22) 



where M e depends on v±, V2, k, I and e, but is independent of R\ and R2 (here we used 
that u = vi + V2 + 2). Since by hypothesis v±, V2, k and / are not all equal to zero, the 
exponent of the integrand in (4.22) is nonnegative and we obtain 



R1R2 

On the other hand, by (4.14) and (4.16), we have that 
AA" 



(4.23) 



h < 



(R 1 -k)(R 2 -l) 



f (x/2)-^+ k+1 +^(x/2)-^ +l+ ^x 2 ( k+l+u -Vdx 

Jl-e 



< 



M 



1 rl 



< 



R1R2 Ji 

eM' 
R1R2 ' 



x V 1 +V 2 +k+l dx 



(4.24) 
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where M' is independent of e, R\ and R2, depending just on v\, v%, k and I. 

By (4.21), (4.23) and (4.24), for R 1 = R + r and R 2 = R the first term in (4.21) is 
majorized by 

eM' M £ 1 / r— 



Extracting the square root we obtain 



Jo 



E 



(v 2 +l,l—V2—l) 

R2-1 



(l-x/2) 



X 



k+l+u-1 



dx 



< 



V7W 



+ 



R ' R 2 ' 

Multiplying the previous inequality by R, we obtain that the quantity whose limit is taken 
in (4.18) is majorized by yeM' + \ / M e /R. This quantity can be made arbitrarily small 
by first choosing e so as to make V eM' arbitrarily small, and then requiring R to be large 
enough to make \ / M e /R arbitrarily small. This proves (4.18). 

A similar argument proves (4.19). This completes the proof of the Lemma. □ 



5. The asymptotics of Mi 

The following result will be needed several times during the remaining part of the paper. 

Lemma 5.1. Let a(t), k(t) and h(t) be complex-valued functions that are real for real t 
and analytic in a domain TcC containing the interval (0, 1] . Let q > be fixed. Then 



[ t Rq h(t) cos[Ra(t) + k(t)]dt 
Jo 



h(l) 



iV<7 2 + (a'(l)) ; 



cos 



Ra(l) + k(l) — arctan 



a'(l) 



+ 0(R~ 2 ). 



(5.1) 



Proof We have 

/ t Rq h(t) cos[Ra(t) + k(t)]dt 
Jo 

' e i(Ra(t)+k(t)) _|_ e -i(Ra(t)+k(t)) 



,Rqlnt 



h(t)dt 



= --{ [° e- R ^ lnt -^e lk ^h(t)dt+ f\- R ^ lnt+ia ^e- ik ^h(t)dt) . 
2 Ui Ji J(5.2) 

The asymptotics for large R of each of the two integrals on the last line of (5.2) can 
be found by the Laplace method as it is described for example in §6, Chapter 4 of [O]. 
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Indeed, consider the first integral. The only requirement of the hypothesis of Theorem 
6.1 of [O, p. 125] that needs to be checked is that the real part of the coefficient of —R in 
the first exponential in the integrand attains its minimum at t = 1. This indeed holds, 
since — hit has its minimum at t = 1. The relevant quantities are easily found to be, in 
the notation of the quoted theorem of [O], A = /i = 1, po = — q — ia'(l), qo = e lk ^h(l) 
and p(l) = — ia(l). By that theorem, the value of the integral is aoe~ Rp ^ /R + 0(R~ 2 ), 
where ao = 9o/(a*Po )• Therefore, we obtain by Theorem 6.1 of [O, p. 125] that 



r° 

J e -R[- q \nt-ia(t)] e ik(t) h ^ dt = 



e i(Ra(l)+k(l)) h ^ 1 



ia'(l) 



R 



+ 0(R~ 2 ). 



(5.3) 



The two integrals on the last line of (5.2) are complex conjugates, so the coefficients of 
their asymptotic expansions are also complex conjugates. We obtain from (5.3) that 



e -R[- q ]nt+ia(t)] e -ik(t) h ^ dt 



However, it is readily checked that 



+ 



e -i(Ra(l)+k(l)) h n\ 1 

q - ia'(l) R 



COs(</7 — 6>), 



+ 0(R- 2 ). 



(5.4) 



q + ib ' q-ib + b 2 

where 9 = arctan(6/g). By (5.2)-(5.5) we obtain the statement of the Lemma 

Lemma 5.2. For fixed r, v\ and V2 we have 

I Vl +i, V2 (R + r, R) 

(_4)^i+«2+i r 1 /a ^\ 1 / 2 



(5.5) 



□ 



2Rn 



x 



X 



X V!+V 2 + 1 



r cos 



1 - 



X 



2/1 



dx + oiR' 1 ). 



(5.6) 



Proof. By Lemma 4.2, it is enough to show that J Vl +i jV2 (R + r, R) has the asymptotics 
given by (5.6). Using (4.17), (4.9) and the fact that sin(cos _1 (l - x/2)/2) = (x/4) 1 / 2 , 
cos(cos _1 (l — x/2)/2) = ((4 — x)/4) 1 / 2 and cos(z — nix) = (— l) n cos(z), one obtains that 



+i,„ 2 ( J R + r, R) 



1 '4-^ 1/2 



7T\/R + r — Vi — 1V-R — V2 Jo \ ■>' 

X COS 



X V1+V2 + 1 cog 



-Rcos -1 1 - 



2) 4 

Converting the product of cosines into a sum, we obtain that 



(R + r) cos" 1 (l-|)- \_ 

(5.7) 



dx. 



+ i,„ 2 (-R + r, R) 



2tx^R + r - v ± - ly/R - v~ 2 
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(J1 + J2), 



(5.8) 



where 



Ji = 



1 '4-^ 1/2 



x 



1 M-^ 1/2 



X 



x vi+v 2 + l cog 



r cos 



IV 

2/. 



(5.9) 



(2 J R + r)cos- 1 (l-|)-^ 



By Lemma 5.1, J2 = 0(1/ R). Therefore, (5.7)-(5.9) imply that the asymptotics of 
J Vl+ i jV2 (R + r, R) is given by the right hand side of (5.6). As noted in the beginning 
of the proof, this completes the proof of the Lemma. □ 



We are now ready to give the asymptotics of M 1 . 
Proposition 5.3. For fixed r, v± andv 2 , we have 



M^R + r, R) 



1 1 r 1 



4 - x 



x 



1/2 



x u 1 COS 



r cos 



1 ( 1_ I)] dx + °( R ~ 3 ^ ( 5 - 10 ) 



where u = v\ + v 2 + 2. 

Proof. Using the bounds (4.14) for the Jacobi polynomials, we obtain from (4.8) that 

m r 1 1 1 



lki(Ri, R2) < 



Jo 



y/RT^VRTn Jo (^^^H (x_^+ 

- / X 

1VR2 Jo 



2+1 1 1 
4 



< 



M 



•X 



l +v 2 +k + l 



dx, 



(5.11) 



where the constants M and M' depend just on vi : v 2 , k and /. Therefore, we obtain that 

1 

R 



I kl (R + r,R) = 



< k < v 1 + 1,0 < I < v 2 . 



(5.12) 



Consider now the representation of Mi given by (4.7). For R\ = R + r and R 2 = R, the 
coefficient of each I k i(Ri,R 2 ) with (k,l) 7^ (vi + l,v 2 ) is 0(R~ 3 ). Thus, by (5.12) we 
obtain from (4.7) that 

M 1 (R + r,R) = (R + r,R)+0 (J^J . (5.13) 

By (3.8), (4.2), (4.3), (4.5) and (4.6) we obtain after simplifications that 

2 

Cv 1 + l,v 2 ~ (_ 4 )„ 1+ „ 2 + i • 

Substituting this value into (5.13) and using Lemma 5.2 we obtain the statement of the 
Proposition. □ 



6. The asymptotics of M a , M c and M ac 

Our analysis of the asymptotics of Mi can be repeated for the remaining double sums 
M a , M c and M ac . We obtain the following result. 
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Proposition 6.1. For fixed r, v± andv 2 , we have 



M a (R + r, R) = - 



1 1 

P* Wo 



ii— i 

X COS 



V COS 



(6.1) 



MAH I r.h>) = -i- / x tt_] cos 
R ^ Jo 



[ x u 1 cos rcos 1 (l - ^) - J cte + o(P 2 ) 



(6.2) 



M ac (P + r, P) 



1 1 

P^ 



A-x 



x 



-1/2 



x u 1 COS 



rcos 



^(l- X -)]dx + o{R^) 



(6.3) 



where u = v\ + i> 2 + 2. 

Proof. By (4.1) we can express the 3P 2 's of (3.9) and (3.10) as 

Ul + l 



3-^2 



-Pi + l,Pi + l,vi + § _ x 
|, vi + 3 '4 



1 



- I ) 



fc/vi + n, 3 



(-Pi + l) fc (Pi + l) fc z 

X ; — — 2-^1 



(Vl + 3)fc 



2M+1 fe=0 

Pi + 1 + k, Pi + 1 + k _ x 
vi + 3 + k '4 



4k 



(6.4) 



and 



3-^2 



-P 2 + 1,P 2 + l,v 2 + I ^ 



;,v 2 + 2 



' 4 



-v 2 



VEH)'(*)(-*i).-. 

~ 2^2 ,_ n V f / Z 



X 



(-P2 + l)/(P 2 + l)/a; 



2/ -"2 ; = 
I 



(v 2 + 2)i 
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>Pi 



P 2 + 1 + /, R 2 + 1 + / . x 
v 2 + 2 + l '4 



(6.5) 



By (4.4), the resulting 2 Pi's are expressed in terms of Jacobi polynomials as 

(Pi 



2 Pi 

and 
2 Pi 



-R 1 + l + k,R 1 + l + k x 
Vl + 3 + k '4 



1 - k - 1)1 (vt + k + 2)! („ 1+fc+2 ,fc-„ 1 -i) / _ x\ 
(Pi+^i + 1)! "~ Rl ~ k ~ 1 V 27 



(6.6) 



-P 2 + 1 + /,P 2 + 1 + / x 
v 2 + 2 + l '4 



_ (R2-I- 1)! (vg + j + 1)! ( V2+ i+i,i- V2 ) f _ x\ 
(R 2 + v 2 )\ ^-z-i V 2)' 



(6.7) 

Substituting the expansions (6.4) and (4.3) and the formulas (6.6) and (4.6) into the 
integral representation (3.11) of M a , we obtain that 

m. = "f t ^t^lM^llk^M^i r u(Bu R2h (6 .8) 



j~ Q R 2 (Pi - k) Vl+k+2 (R 2 -l + l)v 2 +l 
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where c' kl depends only on k, I, v\ and V2, and for < k < v i + 1 and < / < V2, 

I' kl (Ri, R*) := f Q P^El*-*-* (l - f ) P^-i 1 ' 1 -^ (l " |) * k+l+U dx. (6-9) 
Similarly, we get from (3.12) that 

M _ c ki R 2 (-Ri)k(Ri)k (-R2 + iMgg + 1); r „ , p D x , filfV . 
^<c- 2^ Z^ - ^ - 7p — Z~Tn Tr — n hiiR^R?), (0.10) 

where depends only on A;, I, v\ and ^2, and for < k < v\ + 1 and < / < V2, 

l' k \( Rl , R 2 ) := jf 1 p(-+^^--- 2 ) (l _ |) Pgw™ (l - |) x*+^cfa. (6.11) 
An analogous calculation yields from (3.13) that 

II 1 + 1 V 2 



1~t (Ri- k) Vl+k +2 (H2 - i)v 2 +i+i 



{-R 1 + l)fc(i?! + l) fc (-i? 2 + l)l(^2 + 1) 

c kl si2Sii 

k=0 1=0 

where c' kl depends only on k, Z, v\ and V2, and for < k < v 1 + 1 and < / < V2, 



l' k URi, R2) := f Q pK'*" ,1 " 1) (l " |) P£3- + i 1 ''- Wa) (l - |) * fe+/+U+1 ^. (6.13) 

As seen in the proof of Lemma 5.2 for the case of Iki, the bounds (4.14) imply that for 
fixed r, vi, V2 and fixed < k < v\ + 1 and < / < V2, the integrals I' kl (R + r,R), 
I kl (R + r, R) and I k [(R + r, R) are 0(1/ R). Indeed, the only change from that case is that 

now the parameters of the Jacobi polynomials Pn a ' (x) that occur are slightly changed. 
The key fact needed to prove (5.12) was that the exponent of x in the last integral of 
(5.11) was nonnegative. However, the analogous exponents for the case of I' k l , I kl and 
I' kl are readily seen to be nonnegative as well (by (4.14), this exponent goes down half 
a unit for each unit of increase in the a-parameter of the Jacobi polynomials that occur; 
the a-parameters of the pairs of Jacobi polynomials Pn'^\x) appearing in I' kl , I kl and 
I' kl are increased by (1, 0), (0, 1) and (1, 1), respectively; the increase in the exponent of x 
in (3.11)-(3.13), namely 1, 1, and 2 units, respectively, makes up for the decrease due to 
the change in the ct-parameters) . 

Using this, it can be shown that, just as it was the case for the expansion (4.7) of M 1: the 
asymptotics of the double sums (6.8), (6.10) and (6.12) for r, v±, V2 fixed and R\ = R + r, 
R2 = -R, R — * 00 are given by the contribution of the terms with (k, I) = (v± + 1, V2). 

To work these out, note first that analogs of Lemma 4.2 hold for I' Vl+1 V2 , jV2 and 

J"h-i,v2' with ^«i+i,« 2 > J v'i+i,v 2 and J vt+i,v 2 defined by replacing the Jacobi polynomials 
in the integrands of the /-integrals by their Darboux approximants. Indeed, the only 
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difference from the calculations in the proof of Lemma 4.2 is that now the a-parameters 
of the pairs of Jacobi polynomials Pn (x) appearing in the /-integrals are increased 
by (1,0), (0,1) and (1,1), respectively. However, just as was the case in the previous 
paragraph, the increase in the exponent of x in (6.9), (6.11) and (6.13) (of 1, 1 and 
2 units, respectively) compensates the decrease due to the change in the a-parameters. 
Therefore, arguments parallel to the ones in the proof of Lemma 4.2 lead to analogs of 
the majorizations (4.22) and (4.24) that maintain the key feature of having non-negative 
exponents of x in the integrand, and thus prove the claimed analogs of Lemma 4.2. 

Second, using these analogs of Lemma 4.2, one can easily deduce analogs of Lemma 5.2, 
yielding 

I' Vl+1}V2 (R + r, R) 



2Rn 



x Vl+V2+1 cos 



(r - 1) cos" 1 (l - |JJ dx + oiR' 1 ) 



(6.14) 



I': i+1 , V2 (R + r,R) 

(_4)«l+«2+2 

= 2R^ 



J x Vl+V2+1 cos [(r + 1) cos" 1 (l - I)] dx + oiR' 1 ) 



(6.15) 



l':' 1+hV2 (R + r,R) 

JO 



2Rtv 



x 



-1/2 



X V!+V2+1 cog 



r cos 



1 (l-^jdx + oiR- 1 



(6.16) 



And third, replacing in (3.11)-(3.13) the expansions (4.2), (4.3), (6.4) and (6.5) and for- 
mulas (4.5), (4.6), (6.6) and (6.7), the constants c' vi+l v2 , c" i+l v2 and c"' i+l v2 are found, 
after simplifications, to be 



-V!+1,V 2 ^_4^«l+«2+2 

2 



(_4)"l+«2+2 

2 



1+V2+3 ' 



Substituting these and (6.14)-(6.16) into (6.8), (6.10) and (6.12) we obtain the statements 
(6.1)-(6.3) of the Proposition. □ 



7. The asymptotics of the correlation uj{r,u) 

Substituting the asymptotics of the double sums Mi, M a , M c and M ac given by Propo- 
sitions 5.3 and 6.1 into the formula (3.3), we obtain the following result. 
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Proposition 7.1. For fixed r, v\ andv 2 , we have 



uj b (R + r, Vl ;R, v 2 ) = ^|SiS ac + S a S c \ + o(i?" 1 ), 



where 



Si 

S a 

s c 

Sac 

and u = vi + v 2 + 2. 



1 M-^ 1/2 



x u 1 cos 



r cos 1 



"0 



dx 



/ x u 1 cos r cos 1 ( 1 - - ) - - 
Jo L V 27 2 

1 ( x - 1) - 



x u 1 cos 



4 - x 



X 



r cos 

-1/2 



7T 
2J 



<ix 
<ix 



x u 1 cos 



r cos 



1 i 1 ~ I) 



<ix 



Remark 7.2. By the above result, for fixed r, i>i and v 2 the asymptotics of uJb(R + 
r, ui; i?, v 2 ) as — > oo depends only on the sum v± + v 2 , and not individually on v\ and 
i»2- This is consistent with the expectation that the quadromer correlation at the center 
should depend only on the separation vector (r, it) . 

In the statement of Theorem 1.1, the coordinates of the separation vector (r, it) are 
related by it = qr + c, where q > and c are fixed rational numbers. When q > 0, the 
asymptotics of uJb(R + r, v±; R, v 2 ) as R — > oo can be obtained from Lemma 5.1. To handle 
the case q = we need the following result. 

Lemma 7.3. Let a, k and h be real-valued functions that are analytic in an open interval 
containing (0, 1]. Assume a'(t) > in (0, 1) and lim t ^ + h(t)/a'(t) = 0. Then 



h{t) cos[Ra(t) + k(t)]dt 



Hi) 

Ra'(l) 



cos 



7V 



Ra(l) + k(l)-- +0(R~ Z ). 

£ - 



(7.1) 



Proof. As in the proof of Lemma 5.1, express the integrand in terms of exponentials as 



h(t) cos[Ra(t) + k(t)]dt 



e iRa(t) e ik(t) h ^ dt+ I e -iRa(t) e -ik(t) h ^ dt 
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(7.2) 



Consider the first integral on the right hand side of (7.2). Make the change of variables 
y = 1 — t to obtain 

e -*Ri<fi)6(t)dt, 



l e iRa(t) e ik(t) h ^ dt 

Jo 



(7.3) 



where 7(2/) = — a(l — y) and S(y) = e lfc ^ 1_y ' ) /i(l — y). These functions 7(2/) and 5(y) are 
readily checked to satisfy the conditions in the hypothesis of Theorem 13.2 of [O, p. 102] 
with i replaced by — i throughout (clearly, by complex conjugation, the statement of the 
quoted theorem remains true when i is replaced by — i throughout; we need to apply this 
modified version of the quoted theorem because its hypothesis requires 7(2/)' > 0; compare 
with the beginning of §13.1 of [O]). Since the functions a(t), k(t) and h(t) of (7.1) are 
analytic at t = 1, it follows that 7(1/) and S(y) are analytic at y = 0. Therefore, the 
exponents A and /j of (13.02) [O] are both equal to 1. Thus, Theorem 13.2 of [O] is 
applicable and, since we are assuming lim t ^ + h(t)/a'(t) = 0, it yields 



/ 

Jo 



-iRy(t) 



5(t)dt 



8{0) e 



-iR-y(0) 



7'(0) iR 



+ o(iT i ) 



(7.4) 



(since i is now replaced by — i throughout Theorem 13.2 of [O]). By (7.3), (7.4) and the 
definition of ^(i) and 8{t) we obtain that 



p ik(l)u(i \ p iRa(l) 

e iRa(t) e ik(t) h{t)dt = ftUJe + o(irl) 



a'(l) 



iR 



(7.5) 



The two integrals on the right hand side of (7.2) are complex conjugates, so the coefficients 
of their asymptotic expansions are also complex conjugates. We obtain from (7.5) that 



e -iRa(t) e -ik(t) h{t)dt = _ e tk ^h(l)e-^ + ! 

a (1) iR 

By (7.2)-(7.6) we obtain the statement of the Lemma. □ 



(7.6) 



We are now ready to prove our main result. 
Proof of Theorem 1.1. Let u = qr + c, where q > and g,cG 
integral Si of Proposition 7.1 becomes 



Si = 



A-t 



1/2 



t qr+ 



c-1 



COS 



r cos 



i-i 

2 



are fixed. Then the 



dt. 



For q > 0, we can apply Lemma 5.1 with h(t) = t c_3 / 2 (4 - t) 1 / 2 , a(t) = cos _1 (l - t/2) 
and k(t) =0 to obtain 



Si 



1 V3 
r VQ 2 + V3 



cos 



arctan — - +o(r _1 ), q > 0. (7.7) 

3 qV3j 
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If q = 0, we have u = c and Si becomes 

Si -- 



t 



t ' 1 cos 



r cos 



1 - 



dt. 



This has the form of the integral in (7.1), with h(t) = t c - 1 / 2 (4-t) 1 / 2 ? a (t) = cos _1 (l-*/2) 
and k(t) = 0. Clearly, a'(t) = {At - t 2 )~ 1/2 > in (0, 1), and it is readily checked that 
lim 4 ^ + h(t)/a'(t) = (in checking this we need to use the fact that c = u> 1). Therefore, 
Lemma 7.2 yields 

iSl = ? cos (y"i) +o(r " 1) ' q = - (7 - 8) 

A conceptual way of viewing (7.7) and (7.8) together is to say that (7.7) also holds in 
the limit q — > + . Similar applications of Lemma 5.1 to the integrals S a , S c and S ac of 
Proposition 7.1 yield, for q > 0, that 



Sa = ~ 

r 
1 



S c = 

Sac 



I 1 f T-K 7T 1 \ _-, 

: cos arctan — — + o(r ) 

r V? 2 + 1/3 V3 2 qVs) K ' 

— ^ +o(r~ 1 ) 



m n 

arctan — = i 

3 2 q Jl) 



r VQ 2 + V3 

1 1 /^7T M . _ h 

— — ; cos arctan — — + o(r , 

r VZ^q 2 + 1/3 V3 9> /3/ 



(7.9) 
(7.10) 
(7.11) 



the formulas also holding, by Lemma 7.3, in the limit q — > + . 

By (1.4), Proposition 7.1, (7.7), (7.9)-(7.11) and the fact that the latter four relations 
hold also in the limit q — > + , we obtain 



w(u, r) = ^ l^i^oc + 5 a iS c 



1 



47r 2 r 2 ( 9 2 + l/3) 
3 

47T 2 (3q 2 r 2 + r 2 ) 

+ o(r" 2 ). 



2/ r7r 

cos arctan 



qV3 



. ,. rn 1 

+ sin arctan — = 

3 q^3 



o(r~ 



+ o(r 2 ) 



47r 2 (r 2 + 3w 2 ) 
This proves Theorem 1.1. □ 



Acknowledgments. I would like to thank Jeff Geronimo for useful discussions and for 
pointing out to me Laplace's method for the asymptotics of integrals, and the referee for 
the careful reading of the manuscript and helpful suggestions. 



27 



References 



M. Ciucu, Plane partitions I: A generalization of MacMahon's formula, preprint 
(available at the Los Alamos archive, at http: / / arxiv.org/ps/math. CO/9808017| ). 



M. E. Fisher and J. Stephenson, Statistical mechanics of dimers on a plane lattice. 
II. Dimer correlations and monomers, Phys. Rev. (2) 132 (1963), 1411-1431. 
G. Gasper and M. Rahman, "Basic hypergeometric series," Cambridge University 
Press, Cambridge, 1990. 

I. M. Gessel and X. Viennot, Binomial determinants, paths, and hook length for- 
mulae, Adv. in Math. 58 (1985), 300-321. 

R. E. Hartwig, Monomer pair correlations, J. Mathematical Phys. 7 (1966), 286- 
299. 

F. W. J. Olver, Asymptotics and special functions, Academic Press, New York, 
1974. 

G. Szego, Asymptotische Entwicklungen der Jacobischen Polynome, Schriften Kd- 
nigsberger Gel. Ges., (1933), 35-112; reprinted in R. Askey (ed.), "Collected papers 
of Gabor Szego," vol. 2, 399-478, Birkhauser, Boston, 1982. 

G. Szego, "Orthogonal polynomials," American Mathematical Society, Providence, 
R.I., 1975. 



28 



